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ABSTRACT 

The nature of the recently detected HE and VHE 7-ray emission of Tycho's supernova remnant 
(SNR) is studied. A nonlinear kinetic theory of cosmic ray (CR) acceleration in supernova remnants 
(SNRs) is employed to investigate the properties of Tycho's SNR and their correspondence to the 
existing experimental data, taking into account that the ambient interstellar medium (ISM) is expected 
to be clumpy. It is demonstrated that the overall steep 7-ray spectrum observed can be interpreted 
as the superposition of two spectra produced by the CR proton component in two different ISM 
phases: The first 7-ray component, extending up to about 10 14 eV, originates in the diluted warm 
ISM, whereas the second component, extending up to 100 GeV, comes from numerous dense, small- 
scale clouds embedded in this warm ISM. Given the consistency between acceleration theory and 
the observed properties of the nonthermal emission of Tycho's SNR, a very efficient production of 
nuclear CRs in Tycho's SNR is established. The excess of the GeV 7-ray emission due to the clouds' 
contribution above the level expected in the case of a purely homogeneous ISM, is inevitably expected 
in the case of type la SNe. 

Subject headings: (ISM:)cosmic rays - acceleration of particles - shock waves - supernovae individ- 
ual(Tycho's SNR) - radiation mechanisms:non-thermal - gamma-rays:theory 



1. INTRODUCTION 

The detection of high-energy (HE; 100 MeV < 
E < 100 GeV) and very high energy (VHE;£ > 
100 GeV) 7-ray emission from supernova remnants 
(SNRs) is extremely important, because it provides 
direct evidence for the acceleration of charged parti- 
cles (atomic nuclei and/or electrons) inside SNRs to 
energies that are comparable to those of the gamma 
rays. Based on such detections one can hope to 
eventually confirm the idea that SNRs are indeed the 
main source of nuclear cosmic rays (CRs) up to en- 
ergies of about 10 17 eV in the Galaxy, as widely ex- 
pected dBell fc LuceMlgOOHlPtuskin fc ZirakashviliH2003t 
Berezhko fc VoUdl2007t iVolk et al 1120081 IBerezhko et all 
20091 ). see however iParizot et all (j2006Q . 

However, even when 7-rays from a SNR are detected, 
it is not a simple matter to determine the types of en- 
ergetic particles that produce them. If the detection is 
in TeV emission only, like in the case of SN 1006 up 
to now (jAcero et al.l l2010h . then one has to find evi- 
dence that the emission has a substantial 7r°-decay 7- 
ray component produced by nuclear CRs, or whether it 
is merely inverse Compton (IC) and/or Bremsstrahlung 
emission produced by accelerated electrons. 

At a Galactic latitude of about 15 degre es, SN 1006 lies 
abou t 550 pc above the Gal actic gas disk ( Winkle r et al.l 
I2003t iBocchino et~aTl I2011D . From its dipolar symme- 
try in hard X-rays and TeV gamma rays it appears a 
good approximation to assume an interstellar environ- 
ment that is uniform both in gas density as well as in 
magnetic field direction and strength, even though this 
simplification is still in need of investigation. The dipolar 
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symmetry can be attributed to a corresponding charac- 
teristic of the injection process of suprathermal nuclei 
from the shocked downstream plasma into the accelera- 
tion process: such particles can outrun the shock only 
if the outer SNR shock is quasi-parallel to the upstream 
magnetic field and only then can they participate in the 
diffusive shock a c celeration process. A s argued earlier 
(jVolk et al.M2003t IBerezhko et~aT1l2009h this means that 
only approximately 20 percent of the total shock surface 
will contribute to nuclear CR production from a type la 
SNR and nonlinear shock-modification effects like mag- 
netic field amplification in a type la SN, whereas about 
80 percent of the shock surface correspond to a quasi- 
perpendicular shock without significant injection of nu- 
clei and to suppressed acceleration. A spherically sym- 
metric acceleration model therefore has to be "renormal- 
ized" in the sense that the integrated particle spectrum 
obtained in spherical symmetry needs to be diminished 
by a factor w 0.23 Making this correction, the sum of 
the multi-wavelength observations and their theoretical 
interpretation suggests that nuclear CRs are produced in 
SN 1006 wit h high efficiency , as re quired for the Galactic 
CR sources (jBerezhko et al.ll2012D . 

With significant detections at GeV as well as at TeV 
energies of otherwise simple objects like remnants of type 
la supernovae in the Galactic gas disk, there is still 
the question whether the circumstellar environment is 
uniform in gas density or not, and what role this non- 
uniformity plays for the overall observed 7-ray spectrum 
and morphology. This question is discussed here for 
Tycho's SNR. The remnant is spatially unresolved in 
7-rays but otherwise very well studied. The spatially- 



3 This argument should hold also for the probably significantly 
more complex and smaller-scale inhomogeneous but topologically 
equivalent external magnetic field configuration, expect ed for type 
la SN Rs in general and for Tycho's SNR in particular ijVolk et al.l 
12003? ). 
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inte grated TeV-spectru m, detected by the VERITAS ar- 
ray (jAcciari et alj|20ill) , is quite compatible with a theo- 
retical model for a type la - explosion in a strictly uniform 
Interstellar Medium (ISM), s tudied in detail in a previ- 
ous paper (|Volk et al.1 12001 which is in the following 
referred to as VBK08. However, the recent Fermi Large 
Area Telescope (Fermi LAT) detection of hi gh energy 
7-rays above 400 MeV ([Giordano et al.|[2012h disagrees 
with this si mple nonlinear model, sho wing a significant 
GeV excess. iMorlino fc CaprioH (|20I2t ) (in the following 
referred to as MC12) have attempted to understand this 
result assuming a spectrum of CR protons N oc e~ 7 with 
a spectral index 7 = 2.2. This spectrum is considerably 
steeper than the spectrum predicted by VBK08, that im- 
plies 7«2. To obtain this result, MC12 assumed Bohm 
diffusion for all the accelerated particles. However, as 
will be shown below, such an interpretation contains an 
internal contradiction: Bohm diffusion even at the high- 
est particle energies involved is inconsistent with such a 
steep proton spectrum. 

In this paper a new interpretation of the detected 7- 
ray spectrum of Tycho's SNR will be given. It is based 
on the expectation that the actual interstellar medium 
(ISM) is clumpy instead of be ing purely homogeneous 
(|Fieldll96arWolfire et al.ll2003l) . 

2. HOMOGENEOUS ISM 

The present form of the solutions of t he nonlinear ac- 
celeration equations, considered here (see [Bcrczhkol l2008l 
for a review), assumes spherical symmetry. In this ap- 
proximation it is possible to predict the temporal and 
radial evolution of gas density, pressure, and mass veloc- 
ity, together with that of the energy spectrum, as well 
as the spatial distribution of CR nuclei and electrons, 
including the properties of their non-thermal radiation. 

This theoretical model has been used in detail to in- 
vestigate Tycho's S NR as the remnant of a type la SN 
(|Krause et al. I [20081 ) in a homogeneous ISM, in order to 
compare the results with the existing data (VBK08). It 
was argued that consistency of the standard value of stel- 
lar ejecta mass M e j = 1.4M© and a total hydrodynami- 
cal ex plosion energy E sn = 1.2 x 10 51 erg (Badenes et al. 
|2006|) with the gas dynamics, acceleration theory and 
the existing 7-ray measurements required the source 
distance d to exceed 3.3 kpc in order to be consistent 
with the existing HEGRA upper limit for TeV 7-ray 
emission. The corresponding ambient gas number den- 
sity TVg = p/m (where p is the gas density and m is 
the proton mass) had then to be lower than 0.4 cm -3 . 
On the other hand, the rather low distance estimates 
from independent measurements together with internal 
consistency arguments of the theoretical model made it 
even more likely that the actual 7-ray flux from Tycho 
is "only slightly" below the HEGRA upper limit. The 
strong magnetic field amplification produced by acceler- 
ated CRs implied a mean field strength of « 400 ^G (e.g. 
iVolk et all 120051) and as such implied in addition that 
the 7-ray flux is hadronically dominated. The shock was 
modified with an overall compression ratio cr w 5.2 and a 
subshock compression ratio a s « 3.6; the latter is consis- 
tent with the observed radio index (Reynolds & Ellison 
1992). The differential proton spectral index was 7 « 2. 

The TeV 7-ray emission from Tycho detected by VER- 
ITAS ([Acciari et al.ll20lH ) corresponds very well to the 
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Fig. 1. — Spectral energy distributions of the 7-ray emission 
from Tycho's SNR as functions of 7-ray energy e 7 , calculated for 
a source distance d = 3.8 pc, together with the experimental data 
obtained by Fermi and VERITAS. Dashed and solid lines represent 
the contribution of the warm-phase ISM and the total 7-ray energy 
spectrum that includes the contribution of the clouds (see main 
text), respectively. 

above expectation. As can be seen from Fig[T] a new 
7-ray spectrum calculated within the kinetic nonlinear 
theory (shown by the dashed line) is well consistent with 
the VERITAS measurements. This new calculation was 
performed following the usual procedure as described in 
VBK08. For the proton injection rate rj = 3 x 10~ 4 this is 
still compatible with the above-mentioned shock modifi- 
cation and softening of the observed radio synchrotron 
emission spectrum. The new distance d = 3.8 kpc 
and the corresponding new ambient ISM number density 
N g = 0.25 cm -3 were taken in order to fit the o bserved 
TeV 7-ray emission (see also lAcciari "eTaII[20Tll) . 

However, as mentioned in the Introduction, the 7- 
ray spectrum measur ed by the Fermi LAT at energies 
400 MeV to 100 GeV ([Giordano et al.ll20T2T ) is consider- 
ably (by a factor 2 to 5) above the value predicted by 
the kinetic theory (see FigJT]). This excess of GeV 7-ray 
emission, when compared with the theoretical predic- 
tions, requires a more detailed consideration of this ob- 
ject and its environment, taking into account new phys- 
ical factors which had been hitherto neglected. 

3. STEEP CR SPECTRUM? 

The first possibility to explain the complete 7-ray data 
might be the search for a solution with an essentially 
steeper CR spectrum for all energies above sa 1 GeV, 
compared with the hard-spectrum solution discussed in 
the last section. Such a soft-spectrum solution was pro- 
posed by MC12. These authors also used a nonlinear 
approach, assuming a homogeneus circumstellar environ- 
ment of the SNR. However, their model differs essentially 
from the present one in terms of internal consistency. 

For some aspects of the nonlinear problem, that re- 
quire a consistent determination of the CR spectrum 
together with the time-dependent dynamics of the ex- 
panding spherical shock, MCI 2 substitute a simplified 
semi-analytical approach. This is not a principle is- 
sue. In addition, there are several basic differences be- 
tween the present approach and theirs. First of all, 
MCI 2 do not permit that some part of SN shock should 
be locally quasi-perpendicular, where the injection of 
suprathermal, ionized nuclear gas particles from down- 



Nature of the 7-ray emission of Tycho's SNR 



3 



stream is strongly suppressed. Their magnetic field 
model is instead topologically equivalent to assuming a 
monopolar magnetic field. As a consequence they sug- 
gest equally efficient injection/acceleration of CRs all 
over the shock surface. Secondly, these authors ignore 
any nonadiabatic gas heating within the shock precur- 
sor. Therefore the magnetic field pressure in the re- 
sulting model plays a more important role than the gas 
pressure for the determination of the shock parameters . 
On the other hand, nu merical simulations (Bell 120041 : 
iZirakashvili et all [2008) of the nonresonant instability, 
caused by CRs within the precursor, clearly demonstrate 
significant gas heating that is roughly consis tent with the 
model considered here (e.g. iBer ezhko 20Ql). 

Under the above assumptions of MCI 2 it is indeed for- 
mally possible to obtain a solution with a steep CR spec- 
trum, consistent with the existing measurements. How- 
ever, apart from the discussed differences between the 
two approaches, the scheme of MC12 in addition con- 
tains the assumption of Bohm diffusion for all particles, 
also for those of the highest energies. For the obtained 
steep spectrum this amounts to an internal contradiction 
which renders such a solution inconsistent. 

To make this clear requires a brief analysis of the 
general properties of CR acceleration in spherical blast 
waves: 

The expanding spherical SN shock produces a power 
law momentum spectrum of accelerated CRs which ex- 
tends up to a cutoff mo mentum determined by the ex- 
pression (Bcrezhkol ll996f ) 

K&W) = R*V S /A, (1) 

where k denotes the CR diffusion coefficient, R s and V s 
are the radius and the speed of the shock respectively, 
and A is a numerical factor whose value depends on the 
shock expansion law. For a rough estimate one can use 
the value A = 10. Note that the above expression is also 
roughly valid for a nonspherical shock, for example for 
the shock transmitted into a dense clump. In this case 
half of the mean size of the clump I, i.e. 1/2, plays the 
role of the shock radius R s . 

When CRs are scattered mainly by Alfven waves, as 
assumed by MC12, the diffusion coefficient can be writ- 
ten in the form 

n(p) = KBohm(p)E B /E w (k = l/p B ), (2) 

where 

KBohm = PBV/3 (3) 

is the theoretically minimal diffusion coefficient, the so- 
called Bohm diffusion coefficient, v and pb are the parti- 
cle velocity and gyroradius, respectively, Eb = B 2 /(8n) 
is the energy density of the large-scale magnetic field, 
and E w (k) is the Alfven wave magnetic energy density 
p er un it loga rithm of wave number k. 

iBel ([200l found that a non-resonant streaming in- 
stability will occur within the precursor of a strong, ac- 
celerating shock. The diffusive streaming of accelerated 
CRs is expected to be so strong in this region that a 
purely growing MHD mode appears with a growth rate 
that is, at least in the outer part of the precursor, larger 
than the growth rate of the well-known res onant Alfvenic 
mode. In addition it was demonstrated (|Reville fc Belli 



l2012f ) that CRs form filamentary structures in the pre- 
cursor due to their self-generated magnetic field. This 
filamentation results in the growth of a long-wavelength 
instability. It is expected that due to these instabilities 
the external magnetic field -Bism is amplified within the 
entire precursor structure up to B > -Bism, since the 
main effect is produced by the most energetic CRs with 
momenta p ~ p max which populate the whole precursor 
diffusively during their acceleration. 

Th e Alfven wave excitation in the shock precursor (| Belli 
119781 ) corresponds to an additional unstable mode. It 
leads, in particular, to a high level of resonantly scatter- 
ing waves. Therefore, the overall field amplification will 
be the result of all instabilities operating in the precur- 
sor. 

These mechanisms should lead to the Bohm limit 
k(p) = KBohm(p) for CR diffusion which is achieved when 
the condition E w (k) = Eb is is reached. 

CRs with the momentum spectrum N cx p~ 7 produce 
an Alfven wave spectrum E w oc k J ~ 2 oc p 2 ~ 7 , where k 
and p are approximately connected by the relation k = 
1/Pb- It can be represented in the form 

E W (P) = E m x (p/p max ) 2 " 7 , (4) 

where E m = -E w (p max ). 

The amplified magnetic field B in the case of a steep 
CR spectrum is produced by CRs with energy p ~ m p c 
because it is these CRs that provide the main contri- 
bution to the overall CR energy content; here m p de- 
notes the proton rest mass. Therefore the Bohm limit 
condition E w (p) = Eb can be fulfilled for p ~ m p c, 
whereas for CRs with p pa 10 6 m p c one would have 

E w pb 10 ( - 2 ~^ x6 E B » 0.1B B , for 7 = 2.2. This means 
that the CR diffusion coefficient of such high-energy par- 
ticles would exceed the Bohm limit by roughly one or- 
der of magnitude. As a consequence, this amplified field 
would play almost no role for CRs with high momenta 
p S> m p c, and therefore the maximal CR momentum 
that can be achieved would be considerably lower than 
10 6 m p c. This is in contradiction to the assumption of 
MC12 that k(p) = AtBohm(p) for the entire range of CR 
momenta p < 10 6 m p c considered. 

This argument demonstrates that a considerable in- 
crease of the maximum CR momentum due to magnetic 
field amplification is expected only in the case of a hard 
CR spectrum N oc p~~* with 7 < 2, where the CRs with 
the highest momenta p ~ p max provide the main con- 
tribution to the overall CR energy content. Fortunately, 
such a spectrum is expected to be produced by SN shocks 
for injection rates that considerably exceed those used by 
MC12. 

4. CR AND 7-ray PRODUCTION IN A CLUMPY ISM 

The physics aspect which is not included in the 
present kinetic and ( "renormalized" ) spherically sym- 
metric model is an essential inhomogencity of the am- 
bient ISM on spatial scales that are smaller than the 
SNR radius. This inhomogencity is not the result of the 
progenitor star's evolution towards the final supernova 
explosion, for example in the form of a wind and a corre- 
sponding modification of the circumstellar environment. 
It is rather an inherent nonuniformity of the average ISM 
on account of (i) the interplay between its radiative heat- 
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ing by the diffu se galactic U V field and the radiative cool- 
ing of the gas (|Fieldlll965l ) and (ii) the stochastic agita- 
tion of the ISM by the mechanical e nergy input and gas 
heating from supernova explosions (jCox fc SmithI IT9741 : 
IMcKee fc Ostrikerl [l977h . The first effect is a thermal 
instability and thus a mechanism for small-scale cloud 
for mation in the ISM driven by runaway radiative cool- 
ing (|Field et al.lll969D . Specifically the balance between 
line-emission cooling and gas heating due to the ultra- 
violet background radiation leads to two therm ally sta- 
ble equilibrium ISM phases (|Wolfire et al.ll2003H . One of 
them is the so-called warm interstellar medium with a 
typical gas number density N g i ~ 0.1 cm~ 3 and temper- 
ature T\ 8000 K, the other one a cold neutral medium 
with 7Vg2 ~ 10 cm~ 3 and T 2 ~ 100 K. According to sim- 
ulati ons the scale of dense clouds is typically l r = 0.1 pc 
(e.g. lAudit fc Hennebelldl2?M[Mrtllnou"e et all 12012ft 
0. The second effect is a general compressible turbulence 
of the ISM, at least on scales in excess of ~ 1 pc, driven 
by the Galactic supernova explosions. According to 
MHD-simulations it involves high- and low-temperature 
gas co mponents out of ionization equilibri um on all larger 
scal es (|de Avillez fc Br eitschwcrdt 2005), for a review 
see iBreitschwerdt et al.l (|2012jT While the overall pic- 
ture of the ISM is clearly not simple, the evolution of 
a young SNR like that of Tycho's SN will encounter a 
single realization of the stochastic ensemble of density 
fluctuations. For the energy spectrum of energetic par- 
ticles, accelerated at the blast wave, small-scale high- 
density fluctuations of the ISM play the most conspicu- 
ous role because they produce mainly particles with en- 
ergies far below the cutoff energies for a uniform circum- 
stellar medium. To estimate the spectral changes due 
to upstream density variations in an analytical model, 
the typical ISM is therefore treated here as a generalized 
two-phase medium, composed of a pervasive warm/hot 
ISM (phase I) - called here for brevity "warm" ISM - 
and small-scale dense clouds (phase II), embeded in this 
warm ISM. 

In order to determine the specific properties of the CRs 
and their nonthermal emission in the case of such a gen- 
eralized two-phase ISM the latter is approximated here in 
a simple form, as a uniform warm phase with gas number 
density 7V g i plus an ensemble of small-scale dense clouds 
with gas number density N g 2- The warm diluted ISM 
phase is assumed to have a volume filling factor Fj w 1, 
whereas the clouds occupy a small fraction of space with 
filling factor F2 "C 1. It is in addition assumed that most 
of the gas mass is contained in the warm phase, which 
means that F\N S \ » F 2 N g 2. Then the SN shock prop- 
agates in the two-phase ISM without essential changes 
compared with the case o f a purely homog eneous ISM 
with number density N g i ([Inoue et al.ll2012l) . Therefore 
it produces inside the phase I of the ISM roughly the 
same amount of CRs and nonthermal emission as in the 
case of a homogeneous ISM. Then one has to estimate 
the additional contribution of the clouds in order to de- 
termine the overall spectrum of CRs. 

The large-scale SN blast wave, interacting with each 
single cloud, produces a pair of secondary transmitted 

4 Gravitationally bound molecular clouds are not considered 
here, because they are large-scale massive elements of the ISM. 
If present, they would require an individual treatment. 



and reflected shocks. The reflected shock propagates in 
the warm isobaric ISM already heated by the blast wave. 
Due to this fact its Mach number is quite low and there- 
fore its contribution to the overall CR production can be 
neglected. 

The size i? S 2 = ^c/2 and the speed V a 2 ~ 
{N g i/N g2 ) 1/2 V S ~ lO -1 ^ of the transmitted shock are 
both considerably smaller than the corresponding values 
of the SN blast wave. Therefore, according to Eq.([T|), and 
using A — 10 for both shocks, the maximum momentum 
of CRs produced inside the cloud 



Pmax2 — Rs2V s 2 / {R s V s )p max i 



(5) 



is much smaller than the maximal momentum p m axi of 
the CRs produced by the SN shock in the warm ISM: 
Pmax2 "C Pmaxi- Since the ram pressure pV 2 is expected 
to be the same in both cases, and since the amplified 
magnetic field pressure reaches roughly the same fraction 
of the ram pressure, the magnetic field values are roughly 
the same in these two cases: B2 ~ B\. 

An estimate of the CR spectrum produced by the 
transmitted shock starts from the expression for the CR 
distribution function 



(6) 



which is valid at all momenta p > pi n j up to the cutoff 
momentum p max in the case of an unmodified shock, and 
roughly valid within the momentum range p m j < p < 
10m p c of the subshock in the case of a modified shock. 
Using this expression for the case of the SN blast wave 
and for the transmitted shock, the ratio of the two cor- 
responding distribution functions can be found as 



h 



<]2 



iVg2 

N gl 



l-(92-3)/2 



V 

Pinjl 



(7) 



in pi n j < p < 10ro p c, where also the expre ssion for the 
injection momentum (jBcrczhk o et al.lll996f) 

Pinj « rn p V s . (8) 

was used. 

According to the calculation from section 2, the power 
law spectrum with the index qi ss 4.3, determined by 
the subshock compression ratio, extends up to the CR 
momenta p ~ 10m p c. 

Since the cutoff momentum p maX 2 <C Pmaxi of the CR 
spectrum, produced by the transmitted shock, is much 
lower than the corresponding value for the blast wave, 
Pmaxi ~ 10 6 m p c, one can neglect the modification of the 
transmitted shock. This leads to 92 ~ 4. 

Taking into account that the 7-ray production is pro- 
portional to the gas density, for e 7 > 1 GeV one can 
write the relation between the fluxes of 7-rays produced 
due to the two shocks considered: 

F^ 2 (ej) = aF 2 F 7 i(e 7 ) exp(-e 7 /e 7max2 ). (9) 

Here the factor a is determined by the expression 



Q2 ( N g2 \ 1,5 ^10_c\ a3 



V s ) 



(10) 



and the 7-ray cutoff energy is e 7max2 ~ 0.1p maX 2C since 
on average the energy of the 7-rays resulting from in- 
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elastic proton-proton collisions is about one tenth of the 
proton energy: e 7 ~ O.lpc. 

Substituti ng the values of the SN shock speed V* = 
5000 km/s (IKatsuda et alJ 1201(1 IVolk et al]l2008h . the 
number density for the warm phase of the ISM N g i = 
0.25 cm -3 , as well as suitable fit parameter values for 
the cold phase of the ISM in the form of AT g2 s» 23N g i = 
6 cm -3 and F 2 = 0.005, results in p max2 — 10~ 3 p ma xi = 
10 3 m p c and a w 800. It is noted here that such parame- 
ter values for the cold ISM phase correspond rather well 
to th e results of nume rical modelling of the two-phase 
ISM (jlnoue et al.ll2"012l ). 

Then the flux of 7-rays produced inside the clouds can 
be written in the form 

F 72 (e 7 ) =4 J F 1 7 i(e 7 )exp(-e 7 /100GeV). (11) 

The total 7-ray flux F 7 = F 7 i + F l2l expected from 
Tycho's SNR for a two-phase ISM is shown in FigfTJ One 
can see that it fits the existing data in a satisfactory way. 
Note that the considerable increase (by a factor of 5) of 
the 7-ray emission at energies e 7 < 100 GeV over and 
above the case of a purely homogeneous ISM is due to 
the contribution of clumps which contain only 10% of the 
ISM mass. 

5. DISCUSSION 

Small-scale dense clumps of sizes l c <C 0.1i? s can be 
also produced by the accelerating CRs themselves within 
the pr ecursor as the result of the so-called acoustic insta - 
bility (IDorfi fc Drurvlll98l IDrurvl[l98l IBerezhkdll~986T) . 

Small-scale bright structures of angular size 10" were 
recen tly detected in nonthermal X-rays (|Eriksen et al.l 
l20H . For a source distance d = 3.8 kpc the correspond- 
ing spatial size is / « 0.2 pc, which would be roughly con- 
sistent with the sizes of the expected clumps. However, 
the acceleration in such clouds is not expected to reach 
electron energies in the TeV range which could lead to 
synchrotron X-ray emission. Therefore these small-scale 
X-ray structures can not be considered as an indication 
that dense gas clumps of size / c ~ 0.1 pc indeed exist 
inside Tycho's SNR. Their existence rather derives from 
the general properties of the ISM, as discussed above. 

The contribution of dense gas clumps in different 
kinds of emission can be roughly estimated as follows. 
First, consider the thermal X-ray emission. Since, be- 
sides other factors, the flux of thermal X-ray emission 
Fx oc MgTVg is proportional to the gas density N g and 
the total gas mass M g of the source, we have Fx2 /Fxi = 
(F 2 /F 1 )(N g2 /N gl ) 2 w 2.6. Since the temperature dif- 
ference in the two gas phases is not an essential fac- 
tor for the soft X-ra ys with energies below 2 keV (e.g. 
lHamilton et al.lll983[ ) we conclude that the soft thermal 
X-ray emission should be dominated by the contribution 
of dense gas clumps. In the hard X-ray range above 
2 keV, on the other hand, the luminosity is sensitive to 
the gas temperature, roughly as Fx oc T 21 oc V^- 2 (e.g. 
lHamilton et al.l 119831 ) which make the contribution of 
dense gas clumps relatively small due to their lower tem- 
perature. The expected luminosity of individual clumps 
is considerably higher (by a factor of about 500) com- 
pared with the surrounding diluted gas of the same vol- 
ume. However each instrument sees the remnant in pro- 
jection. Therefore the expected ratio of X-ray fluxes from 



the projection volume V = irp 2 L containing the clump 
to the nearby one of the same size which does not contain 
the clump is r = (V + 500V C )/V, where V c = tt/ 3 /6 is 
the clump volume, p is the cross-section of the volume, 
and L is the line of sight length. A maximal value of this 
ratio r « 1 + 330/ c /L is achieved for p = l c /2. It follows 
from this expression that the contrast of X-ray emissions 
varies from r w 4 for the central part of the remnant 
where L w 10 pc to about r « 25 at the limb region with 
L ks 1.4 pc. Therefore we conclude that these clumps 
could be detected in soft X-rays from the limb r e gions . 
In their study with Chandra lCassam-Chena'i et al.l (2007) 
observe some small contribution of thermal X-rays from 
the regions occupied by the swept-up ambient gas. They 
concluded that it is not clear whether the faint lines or 
other residual emission comes from the ejecta or whether 
they arise from shocked ambient medium. In line with 
the idea of the present paper it is suggested here that 
this emission comes from small-scale clouds in the sur- 
rounding ISM. The difficulty is of course that the X-ray 
emission of Tycho's SNR is dominated by the nonthermal 
X-ray emission. In the context of X-ray emission also the 
observed irregular ities in the blast wa ve position around 
the remnant (e.g. iWarren et al.1 120051 ) could at least be 
partly due to the interaction with ambient clouds. 

Secondly, the contribution of dense gas clumps to the 
synchrotron emission of SNRs is estimated here. From 
their measurements of the variations in the expansion 
paramet er in the radio r ange a nd comparison with X-ray 
features [Revnoso et al.l (|1997j ) suggested the presence of 
ambient clouds, shocked by the blast wave. Since the 
radio synchrotron emission is produced by electrons with 
energies less then 1 GeV, Eq.([71) with p — m p c is used 
in order to estimate the ratio of the total synchrotron 
fluxes originating within the two gas phases: F v2 /F v \ oc 
(F 2 /F 1 )(f 2 /f 1 ) ~ 0.1. This shows that the contribution 
of dense gas clumps to the radio emission is expected to 
be small. The luminosity of individual clumps compared 
with the neighbouring region of the same volume is about 
15. Therefore one should be able to detect the clumps 
with an instrument of corresponding angular resolution 
from the limb region where the contrast is expected to 
be r w 2.9. The synchrotron X-ray energy flux scales 
as Fx oc FK cp N g V 2 oc F. Therefore the expected X- 
ray flux from all the clumps within the SNR, Fx 2 ~ 
(F 2 /Fi)Fxi «5x 10 _3 .Fxi, is small compared with the 
total flux Fxi- 

The dense clumps of size Z c 0.1 pc have an angular 
size of about 6". Structures of such sizes (or even smaller) 
near the outer shock of Tycho' s SNR have been stud- 
ied in optical H a emission (e.g. iGhavamian et all 120001 : 
iLee et afl 120101 ) . It is difficult to conclude whether the 
dense clumps can be detected in optical emission, even 
though the lack of smoothness of the optical filaments 
points in this direction. 

It should also be noted that there are some observa- 
tions which show that Tycho's SNR is expanding into 
an ISM with large scale density gradients and is possibly 
interacting with a dense ambient m edium towards the 
northern direction (|Lee et al.l [2004) . If confirmed this 
makes the overall picture of this SNR even more compli- 
cated. It is nevertheless suggested here that the spheri- 
cally symmetric approach assuming a uniform ISM on a 
large scale is roughly valid. 
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According to the above considerations, the excess of 
7-ray GeV emission above the level expected in the case 
of a purely homogeneous ISM due to the contribution of 
small-scale interstellar clouds is inevitably expected in 
the case of type la SNe situated in a relatively dense ISM 
inside the Galactic disk. It is not clear whether such an 
effect is expected in the case of a SNR situated in a much 
more rarefied ISM, like SN 1006, in the uppermost part 
of the Galactic gas disk. On a speculative basis, similar 
effects may also take plac e in Cassiopeia A , where the 
observed 7-ray spectrum (jAbdo et al.l I20l6h looks very 
similar to that of Tycho's SNR. The difference would 
be that in the case of Cassiopeia A numerous observed 
knots from the supernova ejecta might take over the role 
of pre-existing interstellar clumps. 

6. SUMMARY 

The 7-ray spectrum of Tycho's SNR, consistent with 
the measurements by Fermi and VERITAS, is proposed 
to be the superposition of two spectra: the first part, ex- 



tended up to about 100 TeV, is produced by the SN blast 
wave within the dilute "warm" phase of the ambient ISM, 
whereas the second part, with a cutoff at about 100 GeV, 
originates in dense clouds embedded in this warm ISM. 
The remarkable connection between CR production and 
the physical nature of the Galactic ISM becomes evident 
through the characteristics of the spatially integrated 7- 
ray emission of the SNR sources. 



This work has been supported in part by the Depart- 
ment of Federal Target Programs and Projects (Grant 
8404), by the Russian Foundation for Basic Research 
(grants 10-02-00154 and 11-02-12193) and by the Council 
of the President of the Russian Federation for Support of 
Young Scientists and Leading Scientific Schools (project 
NSh-1741. 2012.2). EGB acknowledges the hospitality of 
the Max-Planck-Institut fur Kernphysik, where part of 
this work was carried out. 



REFERENCES 



Abdo, A. A. et al. for the Fermi collaboration. 2010, ApJ710, L92 
Acciari, V. A. et al. for the VERITAS collaboration. 2011, 

ApJ730, L20 
Acero, F., et al. 2010, A&A, 516, A62 

Audit, E. & Hennebelle, P. 2008, ASP Conf. Series, 385, 73 
Audit, E. & Hennebelle, P. 2010, A&A, 511, A76 
de Avillez M.A. & Breitschwerdt, D. 2005, A&A, 436, 585 
Badenes, C, Borkowski, K.J., Hughes, J.P., et al. 2006, ApJ, 645, 
1373 

Bell, A.R. MNRAS, 182, 147 

Bell, A.R. & Lucek, S.G. 2001, MNRAS, 321, 433 
Bell, A.R. 2004, MNRAS, 353, 550 
Berezhko, E.G. 1986, Sov. Astron. Lett., 12, 352 
Berezhko, E.G. 1996, Astropart. Phys., 5, 367 
Berezhko, E.G. 2008, Adv. Space Res., 41, 429 
Berezhko, E.G., Elshin, V.K., & Ksenofontov, L.T. 1996, JETP, 
82, 1 

Berezhko, E.G. & V61k, H.J. 2007, ApJ, 661, L175 
Berezhko, E.G., Ksenofontov, L.T., & Volk, H.J. 2009, A&A, 505, 
169 

Berezhko E.G., Kse nofontov, L.T., & Volk, H.J. 2012, ApJ(in 

press) (arXiv:1209.2789iT) 
Bocchino, F., Orlando, S., Miceli, M., & Petruk, O. 2011, A&A, 

531, A129 

Breitschwerdt, D., dc Avillez, M.A., Feige, J., et al. 2012, AN, 
333, 486 

Cassam-Chenai, G., Hughes, J. P., Ballet, J. et al. 2007, ApJ, 665, 
315 

Cox, D.P. & Smith, B.W. 1974, ApJ, 189, L105 

Dorfi, E. A. & Drury, L. O. C. 1985, Proc. 19th Int. Cosmic Ray 

Conf. (San Diego), Vol. 3, 115 
Drury, L. O. C. 1984, Adv. Space Res., 4, 185 
Eriksen, K. A., Hughes, J. P., Badenes, C. et al. 2011, ApJ, 728, 

L28 

Field, G. B. 1965, ApJ, 142, 531 



Field, G.B., Goldsmith, D.W., & Habing, H.J. 1969, ApJ, 155, 
L149 

Ghavamian, P., Raymond, J., Hartigan, P., & Blair, W.P. 2000, 
ApJ, 535, 266 

Giordano, F., Naumann-Godo, M., Ballet, J., et al. 2012, ApJ, 
744, L2 

Hamilton, A.J.S., Sarazin, C.L., & Chevalier, R.A. 1983, ApJ, 51, 
115 

Inoue, T., Yamazaki, R., Inutsuka, S., & Fukui, Y. 2012, ApJ, 
744, 71 

Katsuda, S., Petre, R., Hughes, J. P., et al. 2010, ApJ, 709, 1387 
Krause, O., Masaomi, T., Tomonori, U., et al. 2008, Nature, 456, 
617 

Lee, J.-J., Koo, B.-C, & Tatematsu, K. 2004, ApJ, 605, L113 
Lee, J.-J., Raymond, J.C., Park, S., Blair, W.P., et al. 2010, ApJ, 
715, L146 

McKee, C.F. & Ostriker, J.P. 1977, ApJ, 218, 148 

Morlino, G. & Caprioli, D. 2012, A&A, 538, A81; in the following 

denoted as MC12. 
Parizot, E., Marcowith, A., Ballet, J., et al. 2006, A&A, 453, 387 
Ptuskin, V.S. & Zirakashvili, V.N. 2003, A&A, 403, 1 
Reville, B. & Bell, A.R. 2012, MNRAS, 419, 2433 
Reynolds, S.P., & Ellison, D.C. 1992, ApJ, 399, L75 
Reynoso, E.M., Moffett, D.A., Goss, W.M., et al. 1997, ApJ, 491, 

816 

Volk, H.J., Berezhko, E.G., & Ksenofontov, L.T. 2003, A&A, 409, 
563 

Volk, H.J., Berezhko, E.G., & Ksenofontov, L.T. 2005, A&A, 433, 
229 

Volk, H.J., Berezhko, E.G., & Ksenofontov, L.T. 2008, A&A, 483, 

529; in the following denoted as VBK08. 
Warren, J.S., Hughes, J. P., Badenes, C, et al. 2005, ApJ, 634, 376 
Winkler, P.F., Gupta, G., & Long, K.S. 2003, ApJ, 585, 324 
Wolfire, M.G., McKee, C.F., Hollenbach, D., & Tielens, 

A.G.G.M. 2003, ApJ, 587, 278 
Zirakashvili, V.N., Ptuskin, V.S., & Volk, H.J. 2008, ApJ, 678, 

225 



